Heavy-Light Few Fermion Clusters at Unitarity 
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We examine the physics of two, three, and four heavy fermions interacting with a single light 
fermion via short-range interactions. Four-particle bosonic Efimov states have proven important 
experimentally and also been the subject of significant theoretical effort. Similar fermionic systems 
are just now being investigated. We find that with some simple interactions the four- and five- 
particle states collapse to the interaction range at smaller mass ratios than the three-body state, 
and also before larger clusters can collapse. These states and their excitations can be studied in 
cold atom experiments, providing unique insights into the role of few-body systems in many-body 
physics. 
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The few- and many-body physics of heterogeneous 
mixtures of cold Fermi atoms is of great interest both 
experimentally [H 0] and theoretically For many- 
body systems, different mass ratios may allow for 'exotic' 
superfluid phases, or phases of mixtures of different con- 
densates to be observed experimentally. Few-body mix- 
tures of heavy and light fermions also allow us to study 
the role of new classes of fermionic few-body resonances, 
in particular searching for regimes in which four- and 
five-particle states can play a significant role. 

The impact of few-body resonances in cold atomic 
gases has been confirmed recently for bosons, where Efi- 
mov states have been predicted theoretically and found 
experimentally [l^ - flit . These states can have significant 
effects on loss rates in cold-atom experiments, and are of 
course very intriguing in their own right. In the present 
work, we examine potential few-body fermion states and 
their impact. Such states are clearly very different from 
the bosonic Efimov states, but offer a potentially intrigu- 
ing prospect for future studies. 

In the many-body heterogeneous system, the simplest 
phase diagram at unitarity for a mass ratio relevant to 
mixtures of K and Li has been examined in Ref. 17. A 



normal state is found at large polarization, particularly 
where the number of heavy fermions greatly exceed the 
number of light fermions. The resulting energy is very 
small, though, suggesting that intriguing states perhaps 
including multi-particle condensates may arise near this 
regime. 

Theoretical investigations of the few-body problem 
with multiple heavy fermions and a one light fermion 
have already provoked a great deal of interest. Interest- 
ingly, calculations for two heavy and one light fermions 
find regions where collapse is possible at mass ratios sig- 
nificantly below the critical value of 13.6 where collapse 
and associated Efimov states are universal. The few- 
body interactions must be fine-tuned to produce these 
states, and hence the few-body physics is non-universal. 
In this study, Nishida et al., discuss the potential effects 



of these states in many-body experiments, finding the 
potential to conduct examinations of strongly-interacting 
gases of dimers and trimers. Our results indicate that an 
even richer phase structure may be appropriate, as four- 
and five-body resonances can appear in addition to the 
resonances of two heavy and one light particle. In fact, 
these resonances can appear at lower mass ratios than 
the two-heavy one-light (2H1L) system. 

Model Hamiltonian: We investigate this non-universal 
few-body regime by calculating ground states in a model 
system of three-, four- and many-fermions interacting 
with a single light particle via short-range interactions. 
The Hamiltonian considered for one light particle of mass 
m and N heavy fermions of mass M is: 
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where I labels the light particle, h the heavy particles, the 
two-body potential V2 acts only between light and heavy 
fermions, and the three-body potential is a function of 
the separations between the light and the two heavy par- 
ticles. In all calculations the heavy-light effective mass 
/i = Mm. I (M-fm) is used, along with the effective range 
of the two-body interaction Tg/ / to set the scales of en- 
ergy and distance. The two-body interactions are tuned 
to unitarity or infinite scattering length. We use both 
Poschl- Teller and Gaussian two-body potentials: 
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= -2^;f — A2exp[-(Ar)V2], 



where w 0.671 is adjusted to unitarity . We also in- 
clude a three-body interaction with the same range as 



the two-body interaction for each pair: 

V3{r,r') = v° ^A2cosh-2(Ar) cosh-2(Ar'). (3) 

The two-heavy one-hght problem has been investi- 
gated by Nishida, et al.[9], who found a regime for 
8.6 < M/m < 13.6 where the interaction can be fine- 
tuned to produce three-body resonances. In a many- 
body system these states would yield an interacting gas 
of dimers and trimers. In this study they explicitly as- 
sume that collapse of four- and five-body states are not 
favored. However, for specific sets of geometrically sym- 
metric heavy-particle coordinates (equilateral triangles 
and regular tetrahedrons, respectively), Nishida^is'] used 
the Born-Oppenheimer approximation to study the static 
potential between the heavy fermions. Significant addi- 
tional attraction was found for these multi-particle states, 
suggesting it is at possible that the larger clusters may be 
bound at lower mass ratios than the two-heavy one-light 
system. 

Methods: We examine this possibility using Quantum 
Monte Carlo techniques. The calculations employ vari- 
ational states to give a variational upper bound to the 
ground-state energy. We assume a trial wave function 
— 4>L^H, where is a positive definite function 
of the light particle coordinates, and '^h is an anti- 
symmetric state of the heavy particle coordinates. All 
coordinates are measured from the system center-of-mass 
to avoid spurious CM motion. 

The calculations are variational and hence produce an 
upper bound to the ground state of the model Hamil- 
tonian. We believe that these results are likely to be 
accurate in most regimes, as the nodal surfaces = 0) 
are quite simple for these systems. For two, three, or four 
particles it is possible to put all the fermions in relative 
p- waves; larger systems require higher partial waves or 
radial excitations. The simplest possible nodal surfaces 
are, as in the Born-Oppenheimer approximation, inde- 
pendent of the position of the light particle. They are 
simply the projected length, projected area, and volume 
of the line, triangle, and tetrahedron connecting the two, 
three, and four heavy fermions, respectively: 

= (ri-r2))-z. (4) 

= r3a2 X (ri - V2) ■ z. 

= r4,123 • (l"3,12 X (i"i - r2)), 

where r^j^ is used to indicate the relative coordinate of 
particle i from the center of mass of the pair or triplet jk. 
These ground states have angular momentum L = 1, 1, 0, 
respectively, for N = 2,3,4. Slightly lower energies are 
obtained in the calculations reported here by putting the 
heavy particles in s- and p-wave orbitals measured from 
the system CM, yielding additional variational freedom. 

Results: The results for the ground-state energies of 
the two-, three-, and four- heavy, one light systems 
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FIG. 1: (color online) Binding Energies versus mass ratio 
M/m for 4H1L, 3H1L, and 2H1L systems with attractive two- 
body potentials. Energies are given in terms of Eq, energies 
less than zero correspond to systems that will collapse as the 
range of the interaction goes to zero. Filled symbols corre- 
spond to the Poschl- Teller two-body potential and open sym- 
bols to the Gaussian potential. 





M/m = 13 


M/m = 14 


2H1L 
3H1L 
4H1L 


-0.001(1) 
-0.088(1) 
-0.340(1) 


-0.0057(6) 
-0.135(1) 
-0.454(1) 



TABLE I: Binding energies for the 3H1L and 4H1L system at 
mass ratios near threshold for the 2H1L system. Energies are 
given in terms of Eq for the Poschl- Teller potential. 

are shown in Figure [T] for these specific two-body in- 
teractions. All energies are plotted in units of Eq = 
fi^/(2/ir^y^), where fe// is the effective range of the in- 
teraction. For these two-body interactions the regimes 
where E <0 are shown, nevertheless they would collapse 
to large negative energies and zero radius as the two- 
body effective range is reduced. The few-particle system 
will collapse at this point independent of whether the 
particles are confined in a trap or not. 

For these interactions we find that the four heavy one 
light system collapses at a much smaller mass ratio than 
the three heavy one light which, in turn, collapses be- 
fore the two heavy one light. The additional attraction 
obtained from the light particle orbiting multiple heavy 
fermions is sufficient to bind them quite deeply. The in- 
crease in binding with mass ratio is very rapid for the 
3H1L and particularly for the 4H1L system. For these 
simple interactions, the binding of the 3H1L and 4H1L 
clusters are very large at the mass ratio where the 2II1L 
system is near threshold. Results for the Poschl- Teller 
potential are given in Table HI 

It could be possible that adding additional fermions 
would bind the system at even smaller mass ratios. To 
test this assertion we have performed calculations of 
the 9H1L and lOHlL systems using the same Quantum 
Monte Carlo techniques. These particle numbers are sug- 





M/m 


VMC L 


VMC H 


QMC L 


QMC H 


2H1L 


14.5 


0.78 


0.18 


1.14 


0.30 


3H1L 


11.75 


0.30 


0.18 


0.42 


0.24 


4H1L 


10.25 


0.30 


0.24 


0.30 


0.24 



TABLE II: RMS radii for the heavy and light particles for 
the clusters. VMC and extrapolated QMC result are given, 
results are for the Poschl- Teller potential and all radii are 
given in units of r^ff 



gested by filling single-particle states in the d (L=2) wave 
and a radial excitation into the 2s state. Adding addi- 
tional fermions does not seem to reduce the minimum 
mass ratio further, presumably the additional energy cost 
of putting heavy fermions in higher angular-momentum 
or higher-energy radial states is too large to bind the 
overall system. There exists a mass ratio where a sin- 
gle light particle could bind an infinite number of heavy 
particles, however in our calculations to date this ratio 
appears significantly larger than the ratio where smaller 
systems can collapse. 

We have calculated the root-mean-square radius for the 
heavy and light particles for the clusters near threshhold. 
The results for a few cases are summarized in Table |lTl 
For the 2H1L system we find that the light-particle rms 
radius at very weak binding is about 3-4 times the heavy 
particle radius, consistent with simple expectations. The 
difference in radii decreases as we go to the larger sys- 
tems and larger binding, though. In the table we report 
VMC results for the trial wave function ^>t and the QMC 
results. The large extrapolations make very accurate re- 
sults difficult to extract, but the trend is clear. 

We've also examined the single-particle densities of the 
clusters. Because of the nodal structures defined above, 
the heavy particles have very small density near the clus- 
ter's center-of-mass. In the case of the simple trial func- 
tions above, the density is exactly zero at the origin. The 
light particle, in contrast, has a maximum density at the 
origin, preferring to occupy a space in the middle of the 
heavy fermions. 

It is clear from our results that different two-body po- 
tentials with the same effective range give different bind- 
ing energies. In fact the mass ratio where the curves 
intercept E = changes with different two-body interac- 
tions. To further quantify these results, we've also stud- 
ied these systems including a repulsive three-body in- 
teraction V3. The three-body interaction also shifts the 
zero binding energy intercept. These results underscore 
the non-universal nature of this regime. As an exam- 
ple, in Figure [2] we compare results with- and without a 
three-body interaction near threshhold for the 3II1L and 
4II1L systems. Even relatively small amplitudes of the 
three-body interaction shift the curves significantly, as is 
to be expected for a non-universal system. The 2II1L 
systems behaves similarly, with a significant impact of 
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FIG. 2; (color online) The effect of three-body interactions 
on the binding energies of 3H1L. The three-body interaction 
can shift the mass threshhold for collapse significantly. 



the three-body force. All systems can become unbound 
for a repulsive three-body interaction with a strength V3 
of order unity. 

Future Possibilities: It would be interesting to know 
if a single three-body interaction parameter can simul- 
taneously describe the binding of the 2H1L, 3H1L, and 
4H1L systems. This has been found to be the case in 
few-body bosonic systems, or in fcrmionic systems where 
the number of spin degrees of freedom is larger than the 
number of particles. In particular this relation has been 
used to explain the relative binding of the deuteron, 311, 
and 4He nuclei. Of course the interaction here cannot be 
a pure contact interaction but would require momentum- 
dependent terms for the two heavy fermions, which must 
be in an odd relative partial wave, coupled to one boson. 

We are presently investigating the possibility of bound 
excited states in these systems in certain ratios of mass 
regimes, perhaps related by simple scaling laws. The 
large binding of the heavier systems relative to the 
smaller gives some reason to believe this may be pos- 
sible. One can imagine states similar to those found for 
bosons where one heavy particle is in an excited loosely 
bound state. It would have to be in an s-wave or other 
low partial wave to have a significant attraction to the 
remaining N-1 particle cluster. Finally, it could be ex- 
tremely interesting to study the systems described here 
in the many-body context. Optical lattices can be used 
to set the effective mass ratio as described in Q. One 
could in principle set the mass ratio and the relative pop- 
ulation of the heavy and light species to obtain strongly- 
interacting gases of 4H1L and 3H1L clusters. Many other 
possibilities can also be imagined. 
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